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The relaxed geometries and electronic properties of the 
hydrogenated phases of the Si(lll)-7x7 surface are studied 
using first-principles molecular dynamics. A monohydride 
phase, with one H per dangling bond adsorbed on the bare 
surface is found to be energetically favorable. Another phase 
where 43 hydrogens saturate the dangling bonds created by 
the removal of the adatoms from the clean surface is found 
to be nearly equivalent energetically. Experimental STM and 
differential reflectance characteristics of the hydrogenated sur- 
faces agree well with the calculated features. 



The identification of the fascinating 7x7 reconstruc- 
tion on the (111) surface of sihcon has motivated con- 
tinued interest in its structural and electronic proper- 
ties over the years. The Si(lll)-7x7 surface is perhaps 
the system that best exemplifies the important interplay 
between structure and electronic properties in surface 
physics. The work of Binnig et al. using scanning tun- 
neling microscopy (STM) led Takayanagi et al. to 
propose their dimer-adatom-stacking fault (DAS) model 
whose structure is consistent with many experiments . 
The geometric and electronic properties of the 7x7 sur- 
face are well known through a combination of experi- 
ments j|-0j2l|] and realistic calculations — especially the 
first-principles investigations of Stich et al. |Q and Brom- 
mcr et al. [|6[ using large supercells in state-of-the-art ap- 
proaches. 

The structure of the hydrogenated phases of this im- 
portant surface are not as well understood. Studies of 
this surface reaction and associated reconstruction are of 
technological relevance in the chemisorption processes of 
H, NH2, Q and hydrocarbons 0. In addition, the hydro- 
genated surface can be used as reference in clarifying the 
optical response of the bare 7x7 substrate, and is then 
important to have a better understanding of its struc- 
tural and dynamical properties [^,^. From a computa- 
tional standpoint, Si(lll)-7x7:H surfaces are challenging 
due to the large size unit cells. 

Experimental evidence suggests that adsorption of hy- 
drogen on Si(lll)-7x7 occurs in a two-step process 
[p|-pl|. The first step likely occurs when only the 19 
dangling bonds (DBs) per unit cell are saturated, corre- 
sponding to those of the 12 adatoms (ADs), 6 rest atoms 
(RAs) and corner hole atom (CH) of the DAS model, to 



form Si(lll)-7x7:19H. The second step occurs when the 
adatom layer is missing and hydrogen saturates the cor- 
responding 43 DBs to form Si(lll)-7x7:43H These 
findings have resulted in several unresolved issues, such as 
whether the DAS structure remains unchanged when the 
surface is terminated by H atoms, which of the two hy- 
drogenated phases is more stable, and what are the main 
electronic properties of the Si(lll)-7x 7 surface after pas- 
sivation with hydrogen. To date, the efforts to address 
these issues theoretically have been limited to different 
reconstructions of this surface and first-principles 

studies of the surface upon H and NH2 adsorption (using 
a cluster model for the substrate) |lj]. An alternative 
and more reliable approach is to perform molecular dy- 
namics supercell calculations that realistically model the 
surface using the full unit cell of the hydrogenated phase, 
combined with the accuracy of an ah initio treatment of 
the interatomic interactions. 

Here, we present the first ah initio studies to charac- 
terize the Si(lll)-7x7:19H and :43H systems using large 
supercell geometries approaching 400 atoms. Applica- 
tion of this method to elucidate the structural and elec- 
tronic properties of the clean 7x7 surfaces yields results 
in good agreement with previous work. We find that H 
adsorption on the clean surfaces leads preferentially to a 
monohydride structure (Si(lll)-7x7:19H), and is accom- 
panied by considerable surface relaxation, in agreement 
with experiments. The other hydrogenated phase identi- 
fied, the Si(lll)-7x7:43H surface, is found to approach 
energetically the :19H phase, explaining the possible co- 
existence reported in experiments. Moreover, we find 
that the electronic density of states of these H-passivated 
surfaces exhibit no surface states in the gap, while still 
showing strong contributions from back-bonding states 
at low H-coverage, in agreement with STM and differen- 
tial reflectance measurements [9|JlO|]. 

Our computational scheme is molecular dynamics sim- 
ulations based on density functional theory in the lo- 
cal density approximation, using the Harris functional, 
Hamann-Schliiter-Chiang pseudopotentials [ p^ , and a 
minimal s-p basis representation. [A comprehensive dis- 
cussion of this scheme, approximations used, and a de- 
scription of various tests can be found in Ref. 14.] In 
silicon, the 7r-bonded (2x1) reconstruction of the (111) 
surface, the 2x1, p(2x2) and c(4x2) phases of the (100) 
surfaces, as well as the 5x5 DAS structure of Si(lll) 
have all been studied by this technique |lH] . 
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As a test, we started applying this approach to investi- 
gate the structural properties of clean Si(lll)-7x7. The 
substrate is represented by a supercell with ten Si-layers 
plus the adatoms on the slab surface (the initial coordi- 
nates of adatoms and the three uppermost Si layers were 
those of Ref. 6). The model system contains all the fea- 
tures of the DAS configuration (see Fig. |l|a). Periodic 
boundary conditions are imposed parallel to the surface, 
and the bottom Si layer is terminated by H atoms to pas- 
sivate dangling bonds. The equilibrium structure of the 
543 atom slab is then obtained with a dynamical quench- 
ing minimization technique. All the atoms, except for the 
bottom Si layer, were allowed to relax until the forces on 
them converged to < 0.1 eV/A. Only the F point of the 
small Brillouin zone of the supercell is used to sample 
the electronic states. The calculated bondlengths of ADs 
with the first layer atoms are found to be stretched from 
3 to 6% with respect to a typical Si-Si bond length (2.35 
A). The average distance between the AD and the sec- 
ond layer atoms directly beneath them is about 2.44 A, 
except for two adatoms on the faulted region of the unit 
cell (shown as 1 and 2 in Fig. pi), which have a distance 
8% larger than the typical value. The bondlengths be- 
tween RAs and second layer atoms are stretched by 4%, 
as they are 0.5 A away from the first layer. The relative 
heights of AD 3 to 6 are 0.069, 0.055, 0.0 and 0.025 A, 
respectively. These values are in qualitative agreement 
with those from FEED, 0.12, 0.08, 0.0 and 0.04 A [0. 
The second and third layer atoms directly below ADs 
move ~0.5 A towards the bulk, and their bondlengths 
are ~ 2% shorter than the bulk Si-Si bondlength. Com- 
parison of our relaxed atomic positions with those from 
Car-Parrinello work show nearly identical geometry 
p^t . We adopted these equilibrated positions for the ge- 
ometry of the clean surface, as it is in general agreement 
with values reported in previous work ||q,[l6|- |l8| , and then 
proceed to saturate dangling bonds. 

The Si(lll)-7x7:19H is modeled by terminating aU the 
silicon DBs of the bare surface with H atoms (see Fig. 
|l|b) , and all atoms — except for the bottom Si layer — of 
the 366-atom supercell are allowed to relax (for computa- 
tional speed, we take only six Si layers here, plus adatoms 
and chemisorbed hydrogens) [|l9|. The positions of the 
surface Si atoms change when the 19 DBs are terminated 
by H atoms, and the heights of ADs become nearly simi- 
lar. Relative to the clean surface, adatoms and first layer 
atoms move away from the bulk layer, while RAs, CH and 
ADs 1 and 2 move towards the slab, leaving a nearly flat 
top layer jl^ . The ADs BB lengths with first layer atoms 
are slightly shorter (~ 1%) than in the clean substrate, 
while those of the RA with second layer atoms are com- 
parable to the bare surface. The calculated Si-H bond 
distances are 1.54 A, nearly 8% larger than the sum of the 
covalent radii of Si and H (1.43 A), and comparable to the 
value obtained from other calculations (1.502-1.520) [ p2[ . 
The dimer atoms and second and third layer atoms di- 



rectly below ADs remain in their clean-surface positions. 
Our calculation predicts an energy of —0.27 eV/ (unit cell 
area) for the hydrogen-covered surface, relative to clean 
substrate plus gas-phase H2 molecules. This indicates 
that the formation of the monohydride phase on Si(lll)- 
7x7 is stable, in agreement with experiments ||9|- pl|p0[ | . 

The Si(lll)-7x7:43H surface has been suggested to co- 
exist with :19H, based on STM observations ||l^. The 
first two layers of the surface atoms in :43H (including 
surface H) , both on the faulted and unfaulted half of the 
unit cell, are found to resemble the hydrogenated bulk 
terminated Si(lll)-lxl surface. As in the :19H mono- 
hydride phase, dimers on the second layer and a corner 
hole atom are also present (see Fig. |l|c). Although AD 
removal is found to have a considerable activation bar- 
rier, it is easily compensated by the formation of Si-H 
bonds [0. This surface is modeled by removing the 
ADs from the clean slab model and passivating the re- 
sulting 43 DBs with H, allowing this 378-atom supercell 
to relax. The equilibrated structure shows Si-H bonds 
directed normal to the surface and bondlengths compa- 
rable to the :19H phase. The uppermost Si atoms are 
three-fold coordinated to the atoms directly below, with 
bondlengths stretched by 1 to 2% relative to the nor- 
mal Si-Si bond, comparing well with FDA results of the 
hydrogenated 1x1 surface ||l^. The RAs and CH move 
towards the bulk by about 0.5 A, respect to the clean sur- 
face, while the second and third layer atoms, originally 
directly below ADs, move away from the slab and reach 
their typical bulk positions. The Si(lll)-7x7:43H phase 
is foimd to be merely 0.12 eV/(Si-H bond) less stable 
than the :19H monohydride phase, indicating that these 
two structures are almost energetically equivalent. This 
result suggests the energetic may indeed allow these two 
phases to coexist, as reported in Ref. 10. 

The calculated electronic density of states (EDOS) also 
contains important information on these systems. The 
clean Si(lll)-7x7 surface contains midgap states which 
we can clearly identify with surface dangling bonds (see 
Fig. 1). These EDOS features, lying above 0.19 eV and 
from to —0.29 eV, arise from the empty and filled AD 
DBs, respectively. [Energies are referred to the Fermi 
level, which is 0.64 eV above the valence band edge 
in close agreement with the experimental value of 0.6 eV 
in Ref. 3.] The dangling bond states are localized at 
various sites with a predominant pz character. Between 
—0.88 and —0.98 eV, we find states resonant with the 
valence band which involve rest atom DBs, while AD 
backbonds lie between —1.38 and —1.44 eV. Corner hole 
states appear at —1.08 and —1.24 eV. Our results are 
in good agreement with several experiments that have 
reported occupied AD and RA dangling bonds at —0.25, 
and —0.85; AD backbonds and corner hole states at —1.8 
eV; and AD back-antibonds at 0.55 eV 

Figures ^-c exhibit simulated STM images accord- 
ing to the scheme of Tersoff and Hamann pa], with 
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the electronic states obtained in our calculation. The 
panels display constant charge density plots (p ~ 0.01 
eV/A^) showing surface states at various energies below 
the Fermi level. Over a wide range of low bias values, 
for energies close to the Fermi level, we see a clear asym- 
metry between the faulted and unfaulted portions of the 
unit cell (not shown here), in agreement with low-bias 
STM measurements |2l]. Figure |^ exhibits states be- 
tween —0.88 and —1.08 eV, corresponding to DBs of RA 
and CH, while Fig. ^b exhibits BB states lying between 
— 1.38 and —1.40 eV, arising from and Py orbitals of 
the AD bonded to Pz orbitals of the first layer atoms. 
The general features of our calculated STM images for 
various RA/CH and AD are in excellent agreement with 
those obtained by Hammers et al. . 

In marked contrast to the clean surface, the EDOS of 
Si(lll)-7x7:19H has no midgap states, as shown in Fig. 
^. This is intuitively expected, as the DBs responsible 
for these states in the bare surface are passivated by H 
atoms. This finding corroborates EDOS measurements in 
Ref. 10 and 20 where features due to DBs on the bare sur- 
face disappear as the substrate is dosed with hydrogen. 
[Also, theoretical STM images of states near £"„ yield fea- 
tureless isodensity plots.] Moreover, we find that states 
lying between —0.74 eV and £"„ are extended and show 
no contribution from the surface H and ADs states. No- 
tice, however, that backbond states associated with the 
surface adatoms are still present, and lie between 1.07 
and 1.21 eV below the Fermi level, and yield the STM 
image shown in Fig. |^c. This finding supports interpre- 
tation of recent differential reflectance experiments 
as discussed below 

The EDOS of the Si(lll)-7x7:43H also contains no 
midgap states (see Fig. A clear reduction of the den- 
sity of states in the region below and in the con- 
duction band is produced by the disappearance of ADs 
states. [Notice in particular the absence of a sharp fea- 
ture at ^ —1.6 eV.] Our results agree quite well with the 
local tunneling spectroscopy experiment of Mortensen et 
al. ijl^ , who find no high density of states near the Fermi 
level for the :43H-coverage. Only extended states due to 
surface and bulk Si were found between —0.8 eV and the 
valence band edge Ey , as all the dangling bonds are fully 
saturated. 

Our findings also provide strong support for the inter- 
pretation of differential reflectance (DR) measurements 
by Roy and Borensztein [^, who employ this powerful 
technique to explore the evolution of the surface elec- 
tronic states under various chemisorption processes. In 
particular, at zero and low H-coverage, a strong DR fea- 
ture at 1.6 eV is identified as coming from transitions 
involving the backbond states of adatoms, as supported 
by tight- binding calculations of model structures [^. As 
the H-coverage increases, however, this feature saturates 
and higher energy transitions dominate the DR spectra. 
This behavior is consistent with the disappearance of the 



backbond states as the higher coverage is reached and 
the :43H phase forms. Notice that the saturation of the 
1.6 eV feature is accompanied in experiments by a strong 
reconstruction of the surface, as probed by LEED, and 
as expected from the relaxation pattern shown in Fig. 
|l|c 0. Finally, we should also notice that an empirical- 
parameter calculation of the optical response of this sur- 
face using our fully relaxed atomic positions yields ex- 
cellent agreement with the DR experiments, as will be 
discussed elsewhere [ p3[ . 

In summary, we obtained significant geometrical relax- 
ation arising from H- adsorption on the clean 7x7 surface, 
accompanied by strong shifts of the character and ener- 
gies of the electronic states. These changes are found 
to be in excellent agreement with STM, differential re- 
flectance, and LEED experiments, and provide the basis 
for detailed comparisons with new experimental measure- 
ments. 
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FIG. 1. Top view of (a) the DAS model for the Si(lll)-7x7 
bare surface (big gray, black and small gray circles represent 
adatoms (AD), rest atoms (RA) and corner hole atom (CH), 
and first layer atoms, respectively); (b) the :19H coverage 
(small white circles denotes H atoms); and (c) the :43H cov- 
erage surface. The faulted portion of the outlined unit cell is 
on the left. Unit cell shown by dashed lines. 

FIG. 2. Electronic density of states for the clean 
Si(lll)-7x7 surface, as well as for the two hydrogenated sur- 
faces studied with 19 and 43 H-atoms per cell. All energies re- 
ferred to the Fermi level; level broadening used 0.1 eV; curves 
displaced vertically for clarity. 

FIG. 3. Simulated STM images in the Si(lll)-7x7 system, 
(a) RA/CH levels (at 1 eV in Fig. 2), and (b) AD back- 
bond states (~ —1.4 eV) for the clean surface, (c) AD back- 
bond states in the monohydride :19H phase (~ —1.2 eV in 
Fig. 2). The faulted portion of the outlined unit cell is on the 
right. 
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Electronic Density of States 




